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A.  Technical  Problem 

In  recent  years  it  has  been  noted  that  small  additions  of  indium  (on 
the  order  of  1  atomic  %)  to  gallium  arsenide  grown  from  the  melt  greatly 
reduces  the  dislocation  density  to  less  than  10  cm  .  The  mechanism  by 
which  this  phenomenon  occurs  has  been  speculated  to  be  solid  solution 
strengthening,  in  which  In  together  with  its  four  nearest  neighbor  As  atoms 
acts  as  the  hardening  unit.  In  order  to  verify  this  strengthening 
mechanism,  a  program  has  been  undertaken  to  determine  the  strength  of 
Ga^_.cInJtAs  alloys  as  a  function  of  temperature  and  composition.  In  the 
first  year  of  the  study,  Vickers  hardness  measurements  as  a  function  of 
temperature  were  used  as  a  preliminary  assessment  of  solid  solution 
strengthening,  since  microhardness  is  a  good  first  cut  to  examine  if 
materials  exhibit  a  plateau  stress  associated  with  solid  solution 
strengthening.  These  experiments  required  construction  of  a  hardness 
tester  described  in  an  earlier  report.  This  second  year  of  study, 
consisting  of  compression  tests  and  transmission  electron  microscopy,  will 
be  described  herein.  In  addition  to  work  on  III-V  system,  it  is  of 
interest  to  pursue  similar  studies  in  II -VI  compounds  which  potentially 
demonstrate  enhanced  mechanical  properties  by  the  addition  of  isoelectronic 
elements.  Preliminary  studies  on  the  II-VI  compound  Cd^_xMnJ{Te  will  also 
be  described. 

B.  General  Methodology 

The  experiments  used  to  study  solid  solution  strengthening  in  the 

second  phase  of  this  work  are  uniaxial  compression  tests.  Compression 

20  _  3 

testing  of  undoped  GaAs  and  GaAs  containing  1-2x10  atoms -cm  of  In  were 
carried  out  at  temperatures  of  500,  700,  900,  1000  and  1100°C  at  a  strain 
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rate  of  10‘^s"^  in  the  [C01]  orientation  favoring  the  operation  of  multiple 
slip  systems  and  in  the  [123]  orientation  favoring  a  single  slip  system. 

Compression  testing  as  a  function  of  temperature  in  single  slip  and 
multiple  slip  orientations  were  used  to  determine  the  critical  resolved 
shear  stress  (CRSS)  directly  from  the  stress  strain  curve.  In  these  tests, 
the  specimen  surface  is  protected  in  a  bath  of  to  prevent  As 

volatilization  and  this  allows  tests  to  be  performed  at  temperatures  near 
the  melting  point.  These  tests  are  more  difficult  to  corduct  than  hardness 
tests  but  supersede  the  problem  of  surface  degradation  present  in  the  hot 
hardness  tests  at  elevated  temperatures. 

Additional  work  on  strain  rate  dependence  has  been  carried  out  on  a 
limited  basis.  Stress  relaxation  measurements  were  made  following 
compression  tests  in  the  [001]  orientation  at  700,  900  and  1100°C.  Because 
significant  microyielding  is  expected  at  stresses  below  the  yield  point  at 
these  high  temperatures,  it  is  important  to  look  at  the  dislocation 
structure  development  as  a  function  of  stress  level,  m  stly  around  the 
macroscopic  yield  point.  To  start  with,  In-doped  specimens  were  stressed  to 
above  and  below  the  yield  point  at  700°C  at  a  strain  rate  of  10‘4  s  . 

Dislocation  structures  of  undeformed  and  deformed  specimens  have  been 
characterized  by  transmission  electron  microscopy.  Specimens  of  GaAs  and 
Gag  gglnQ  qiAs  f°r  TEM  observations  were  made  by  cutting  thin  slices  in  a 
desired  orientation,  followed  by  mechanical  polishing  and  ion-milling. 

These  specimens  were  examined  in  a  JE0L  JEM  200CX  transmission  electron 
microscope  at  200KeV  using  a  double  tilt  goniometer  stage.  Characterization 
of  dislocation  structures  has  been  performed  by  the  standard  g*b  criterion. 
Electron  microscopy  results  together  with  deformation  data  have  been  used 
to  identify  the  mechanism  of  solid  solution  strengthening. 

Crystals  of  compositions  across  the  CdTe-MnTe  phase  diagram  were 


obtained  from  Professor  Jack  Furdyna  of  University  of  Notre  Dame .  Small 
crystal  wafers  of  Cd^^Mr^Te  (0  <  x  <  0.6)  about  1  mm  thick  were  cut 
parallel  to  the  cleavage  plane  (110),  polished  first  mechanically  down  to 
0.25  micron  diamond  and  then  chemically  in  a  solution  of  l^CroO^-HNO^-l^O. 
The  orientation  was  identified  using  cleaved  edges  and  Vickers  indentations 
were  made  with  one  of  the  diagonals  parallel  to  [100]  direction. 

C.  Technical  Results: 

Compression  Tests 

The  results  of  uniaxial  compression  tests  are  summarized  and  compared 

with  data  from  other  investigations  in  Figures  1  and  2.  Our  results  show 

that  the  critical  resolved  shear  stress  (CRSS)  nearly  doubled  with  In 
20  -  3 

additions  of  1-2x10  atoms  cm  .  Also,  the  CRSS  of  both  undoped  and  In- 
doped  GaAs  showed  a  weak  temperature  dependence,  as  expected  from  the  solid 
solution  strengthening  model.  The  results  also  showed  that  the  onset  of 
dynamic  recovery  process  in  In-doped  GaAs  occurred  at  higher  stress  and 
strain  levels  suggesting  that  the  climb  process  is  more  difficult  in  In- 
doped  alloys  compared  to  the  undoped  alloy.  The  strengthening  observed 
seems  sufficient  to  prevent  macroscopic  yield  during  the  growth  of  large 
diameter  GaAs  crystals  on  the  basis  of  current  models  of  thermal  stresses 
during  crystal  growth.  These  results  and  their  interpretation  are  reported 
in  detail  in  Appendix  A. 

The  strain  rate  versus  flow  stress  curves  that  were  obtained  also 
showed  enhanced  creep  resistance  by  In-doping.  Compression  tests  involving 
abrupt  changes  in  strain  rate  were  performed  at  700°C,  however,  more  te: ts 
at  900  and  1100°C  are  required  to  provide  a  good  understanding  of  the 
strain  rate  dependence  in  this  temperature  range. 


600  700  800  900  1000  1100  1200  1300  1400  1500 

Temperature ,  K 

Figure  2.  A  comparison  of  the  critical  resolved  shear  stress  values  of 

undoped  and  In-doped  GaAs  as  a  function  of  temperature.  Inverted 
open  triangles,  In-doped  and  inverted  solid  triangle,  undoped 
(this  study);  open  triangle.  In-doped  and  solid  triangles, 
undoped  (S.  McQuigan  et  al.,  1986);  solid  squares,  undoped 
(Tabache  et  al.,  1986);  solid  circles  (Swaminathan  et  al.,  1975). 
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Transmission  Electron  Microscopy 

Electron  microscopy  studies  performed  on  undeformed  specimens  revealed 
very  few  dislocations.  Most  of  the  foil  areas,  particularly  in  the  case  of 
In  doped  GaAs  were  found  to  be  dislocation  free.  If  any  dislocations  were 
found,  they  were  generally  straight.  In-doped  and  undoped  GaAs  specimens 
deformed  at  500-700°C  in  the  multiple  slip  orientation,  [001] ,  as  well  as 


n 


in  a  single  slip  orientation,  [123],  revealed  complex  dislocation 
structures.  The  main  dislocation  structures  were  found  to  be  dislocation 
dipoles  (perfect  as  well  as  faulted),  dislocation  loops,  arrays  of  small 
dislocation  loops,  dislocation  tangles,  dislocation  networks  and  straight 
dislocations.  Figure  3  shows  a  typical  electron  micrograph  from  a  500°C 
undoped  GaAs  specimen  deformed  to  Stage  II.  Dislocation  loops  and  faulted 
dipoles  are  clearly  seen.  Figure  U  shows  an  electron  micrograph  of  In- 
doped  GaAs  deformed  at  700°C  to  Stage  II.  Detailed  analysis  of  dipoles  and 
dislocation  loops  were  done  using  dark-field  weak  beam  pictures  for  ±g 
(220)  type  reflections.  The  break  up  of  dipoles  due  to  the  mutual 
attraction  of  the  positive  and  negative  edge  dislocations  of  its  two 
elongated  sides  results  in  a  row  of  loops.  Dislocation  loops  were  found 
to  be  of  the  vacancy  type  with  Burgers  vector  a/3<lll>.  Further  details 
are  found  in  Appendix  B. 


Specimens  in  both  single  and  multiple  slip  orientations  deformed  in 


the  temperature  range  900-1100°C  have  been  taken  to  stage  III  to  IV. 


Therefore,  dislocation  structures  of  these  specimens  were  found  to  be 


fully  recovered.  Dislocation  loops  were  rarely  seen.  The  main  dislocation 


structures  found  in  these  specimens  were  (i)  parallel  sets  of  dislocations 


having  few  interconnections  (ii)  hexagonal  networks  of  dislocations  and 


(iii)  isolated  dislocations.  Figure  5.  an  electron  micrograph  of  900°C 


deformed  In-doped  GaAs,  shows  extended  dislocations  forming  nodes.  From 


Electron  micrograph  of  undoped  GaAs  deformed  at  500°C  to  3% 
strain.  An  array  of  dislocation  loops  (arrow  )  and  faulted 
dipoles  (F)  are  seen  in  the  image.  Zone  axis  —  [001]. 


Electron  micrograph  of  In-doped  GaAs  deformed  at  700°C  t  10% 
strain.  The  dislocation  structures  are  typical  of  deformation 
Stages  II-III.  Zone  axis  -  [111]. 
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the  weak  beam  pictures  taken  in  different  reflections,  the  majority  of  the 
dislocations  in  Figure  5  were  found  to  be  dissociated  into  two  Shockley 
partials .  Figure  6  shows  parallel  set  of  dislocations  having  few 
interconnections.  This  configuration  of  dislocations  is  an  intermediate 
stage  before  recovery.  In  this  configuration  two  out  of  plane  60° 
dislocations  react  to  form  a  third  dislocation  in  a  screw  orientation.  Th 
latter  then  reacts  with  a  dislocation  already  present  in  the  main  slip 
plane  and  finally  forms  a  hexagonal  network  in  a  screw  orientation.  Figur 
7  shows  such  a  hexagonal  network.  This  type  of  dislocation  structure 
indicates  that  climb  process  is  dominant  in  the  recovery  of  dislocated 
structures.  Through  numerous  observations,  recovery  appears  to  be  faster 
in  the  case  of  undoped  GaAs  compared  to  In-dope  GaAs  on  the  basis  of 
dislocation  structures  present  in  the  two  cases.  In  addition,  dislocation 
structures  were  found  to  be  similar  to  those  observed  for  Stage  IV 
deformation  of  Si  and  Ge ,  a  case  where  climb  is  known  to  be  the  dominant 
recovery  mechanism. 

The  stacking  fault  energy  has  been  estimated  by  measuring  the 
separation  of  dissociated  dislocations  (from  weak  beam  micrographs) .  The 
value  of  the  stacking  fault  energy  was  found  to  be  45.3  mJm  .  No  change 
in  the  value  of  stacking  fault  energy  was  found  due  to  In  doping,  implying 
that  there  should  be  little  influence  on  the  cross-slip  probabilitv. 

High  Temperature  Hardness  of  Cd^^.n^Te 

Results  on  the  hardness  of  five  compositions  of  Cd^_xMnxTe  are  shown 
in  Figure  8,  Of  the  five  compositions,  Cdg  g^MnQ  QgTe  showed  numerous 
twins  during  the  optical  microscopy  examination  and  the  measured  hardness 
may  not  be  representative  of  a  single  crystal.  Interpretation  of  these 
results  is  forthcoming. 
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Figure  5.  Electron  micrograph  of  dislocation  structures  in_ In- doped  GaAs 
deformed  at  900°C  for  5.1%  strain.  Zone  axis  -  [111]. 
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Figure  6.  Electron  micrograph  of  the  intermediate  stage  dislocation 

configuration  prior  to  recovery  in  In-doped  GaAs  deformed  at 
1000°C  to  7.2%  strain.  Zone  axis  =  [Oil]. 
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Figure  7.  Electron  micrograph  showing  well - recovered  hexagonal  networks  of 
dislocations  in  In-doped  GaAs  deformed  at  1000°C  to  7.2%  strain. 
Zone  axis  -  [0111. 
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D.  Implications  for  Future  Research 

Both  the  undoped  and  In- doped  GaAs  tested  thus  far  contained  large 
amounts  of  boron  (»  2x10  atoms  cm  )  which  is  also  expected  to  be  a 
potent  isovalent  solid  solution  strengthener .  To  isolate  the  effects  of 
boron  additions,  tests  of  undoped  GaAs  containing  two  different  amounts  of 
boron  are  required.  Undoped  GaAs  prepared  by  a  modified  Bridgman  technique 
with  varying  boron  content  has  now  been  obtained.  Compression  testing  of 
these  crystals  in  the  [123]  orientation  at  different  strain  rates  and 
temperatures  is  in  progress. 

To  better  understand  deformation  during  crystal  growth,  compression 
tests  above  1100°C  will  be  attempted.  This  should  provide  more  significant 
data  for  extrapolation  to  the  melting  point  of  GaAs. 

More  transmission  electron  microscopy  is  planned  for  detailed  analyses 
of  dipoles  and  dislocation  loops.  Faulted  dipoles  will  be  studied  with  the 
combination  of  weak  beam  and  high  resolution  transmission  electron 
microscopy.  From  the  high  resolution  imaging  of  faulted  dipoles  in  In-doped 
and  undoped  GaAs  a  direct  estimate  of  stacking  fault  energy  can  be 
accomplished.  This  value  of  stacking  fault  energy  will  be  compared  to 
the  value  obtained  by  measuring  the  separation  of  isolated  dissociated 
dislocations.  In  addition,  electron  microscopy  of  high  temperature  deformed 
undoped  GaAs  with  varying  boron  contents  is  planned  in  conjunction  with  the 
compression  tests. 


s 
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Further  hardness  and  compression  tests  on  II-VI  compounds  (Cd^_xMnxTe) 
will  continue  following  compression  tests  on  undoped  GaAs  to  examine  solid 
solution  strengthening  in  a  variety  of  semiconductors. 
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Appendix  A 

To  appear  in  J.  Applied  Physics, 
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ABSTRACT 


Compressive  deformation  of  undoped  and  In-doped  GaAs 
single  crystals  has  been  carried  out  in  [00 1]  and  [123] 
orientations  in  the  temperature  range  700  to  1100  °C.  Indium 
additions,  at  levels  of  1-2  x  1020  atoms. cm”3,  results  in  critical 
resolved  shear  stress  (CRSS)  values  that  are  about  twice  as  large 
as  the  undoped  crystals  in  the  temperature  range  of  700  to  1100°C. 
The  CRSS  was  weakly  dependent  on  temperature  in  the  temperature 
range  investigated  as  expected  for  a  model  of  athermal  solid 
solution  hardening.  The  CRSS  value  of  3.3  MPa  for  the  In-doped 
crystal  is  sufficient  to  eliminate  profuse  dislocation  formation 
in  a  75  mm  diameter  crystal  on  the  basis  of  current  theories  for 
the  magnitude  of  the  thermal  stress  experienced  during  growth. 

The  results  also  suggest  that  the  process  of  dislocation  climb  is 


It? 


I*. 


6 


8 


INTRODUCTION 


U* 


v  m 

V. 

''lr 


£ 


s 


VB 

I 


✓ 


Doping  of  GaAs  with  In,  at  a  level  of  about  5  x  1019- 
lxl020/cm3,  in  single  crystals  grown  by  the  liquid  encapsulated 
Czochralski  process  reduces  the  dislocation  density  from  104- 
105/cm2  to  $  102/cm2.1  The  generation  of  dislocations  during  LEC 
growth  in  GaAs  and  other  Ill-V  compounds  is  believed  to  occur 
when  the  thermal  stress  imposed  on  the  cyrstal  during  growth 
exceeds  the  critical  resolved  shear  stress  (CRSS).  Minimization 
of  radial  thermal  gradients  during  growth  and  enhancement  of  the 
inherent  strength  of  the  crystal  by  the  addition  of  isovalent  as 
well  as  Group  IV  or  VI  elements  can  result  in  low  dislocation 
density  GaAs  crystals.3-5  Of  the  possible  dopants,  In  has  been 
found  to  be  very  effective  and  desirable  because  of  its  minimal 
influence  on  the  electrical  behavior  of  GaAs.  Thermal  stress 
calculations  suggest  that  the  maximum  stress  experienced  is  much 
larger  than  the  extrapolated  CRSS  value  of  an  undoped  crystal.6 
The  In  doping,  therefore,  must  result  in  a  large  increase  in  high 
temperature  strength.  It  has  been  suggested  that  hardening  akin 
to  solid  solution  hardening  occurs  in  Ga1_xInxAs  with  an  InAs4 
tetrahedral  cluster  being  the  solute  unit  that  causes 

,  n 

strengthening.  However,  experimental  data  on  high  temperature 
deformation  of  GaAs  and  the  influence  of  dopants  is  scarce  and 
the  role  of  In  is  still  not  well  understood. 

In  a  preliminary  study,  we  first  measured  the  temperature 
dependence  of  the  hardness  of  undoped  and  In-doped  GaAs,  which 
was  reported  earlier.8  In  this  subsequent  study,  we  have 
evaluated  the  flow  stress  in  the  [001]  orientation  for  which 


multiple  equivalent  slip  systems  of  {111}  <110>  type  operate  and 
in  the  [123]  orientation  suitable  for  the  operation  of  only  one 
slip  system,  (111) [101].  The  results  of  the  study  made  in  the 
temperature  range  700-1100°C  are  reported  here,  compared  with 
other  reported  results, 6,9-13  and  discussed  in  terms  of  models 
for  the  high  temperature  deformation  of  diamond  cubic 
materials . 14 • 15 


EXPERIMENTAL  PROCEDURE 


The  75-mm  diameter  single  crystals  of  semi-insulating  GaAs 
and  Ga0<ggln0>01As  used  in  the  present  study  were  obtained  from 
the  Westinghouse  R  &  D  Center.  The  crystals  were  grown  in  [001] 
axial  orientation.  The  undoped  crystal  had  a  dislocation  density 
of  104-105cm-2.  The  indium  doped  crystal  had  an  In  concentration 
of  about  1-2  x  1020atoms.cm“3.  The  boron  concentration  in  both 
doped  and  undoped  crystals  was  about  5  x  1017atoms.cm“3.  The 
dimensions  of  specimens  for  compressive  deformation  in  the  [001] 
orientation  were  5.3  mm  x  5.3  mm  x  10.6  mm  in  size  with  (110'  and 
(110  .<  lateral  faces.  The  dimensions  for  deformation  in  the 
[123]  orientation  were  2.75  mm  x  2.75  mm  x  5.5  mm.  The  lateral 
faces  were  parallel  to  (111)  or  (541)  planes.  The  specimen 
faces  were  mechanically  polished  to  0.3  *im  alumina  followed  by 
chemical  polishing  in  a  1%  bromine-methanol  solution.  The 
specimen  faces  were  parallel  and  orthogonal  to  within  0.5 
degrees . 

The  compression  testing  was  done  on  an  INSTRON  1322 
servohydraul ic  machine.  A  schematic  diagram  of  the  experimental 


set  up  is  shown  in  Figure  1.  Tests  were  performed  at  temperatures 
of  700,  900,  1000  and  1100°C  at  a  strain  rate  of  1  x  10”4s_1  in 
UHP  argon.  The  furnace  used  in  the  study  had  a  low  heat  capacity 
with  a  time  of  heat  up  to  the  test  temperature  of  about  20 
minutes.  For  tests  at  temperatures  of  900°C  and  above  the 
specimen  was  immersed  in  B203  liquid.  The  compression  rams  used 
were  made  of  high  purity  alumina,  and  a  special  fixture  made  of 
alumina  was  used  for  the  alignment  and  positioning  of  the 
specimen.  The  elongation  was  measured  on  the  alumina  rams  by 
means  of  an  extensoraeter  and  was  corrected  for  the  fixture 
compliance,  measured  at  temperature,  in  converting  the  data  to 
strain. 

Stress  relaxation  tests  were  performed  following  the  [001] 
specimen  compression  tests.  At  the  end  of  compressive 
deformation,  the  decrease  in  stress  at  constant  specimen  length 
was  measured  as  a  function  of  time  and  was  converted  to  plastic 
strain  as  a  function  of  time.  The  total  amount  of  strain  during 
stress  relaxation  is  very  small  and  the  total  time  of  testing  is 
only  a  few  minutes. 

RESULTS 

(A)  [001]  Orientation 

Figures  2a  and  2b  show  the  results  of  compression  tests  in 
'001]  orientation  for  undoped  and  In-doped  GaAs.  No  easy  glide, 
stage  I  deformation  is  observed  for  these  crystals,  oriented  for 
multiple  slip.  Also,  only  at  700°C  for  the  In-doped  case  is  there 
an  indication  of  the  yield  point  phenomena  characteristic  of  a 
similar  material6  tested  at  350  to  590°C.  The  tendency  for  less 
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pronounced  yield  point  with  increasing  temperature  is  consistent 
with  results  in  other  studies,6'10  and  the  absence  of  well- 
defined  yield  points  at  high  temperatures  agrees  with  the  work  of 
Djemel  et  al.11  We  note  that  the  slope  below  the  yield  point  at 
700°C  for  the  In-doped  case  is  less  than  the  reported  elastic 
modulus  so  that  microyield  and  plastic  flow  are  already  occurring. 
This  characteristic  has  been  observed  in  other 
measurements,6'10'11  including  those  at  lower  temperatures. 

Indeed,  in-situ  X-ray  topographic  measurements  made  during  tensile 
testing  of  In-doped  and  undoped  GaAs  at  450-700°C  demonstrate 
directly  that  dislocation  formation  and  propagation  occur  at 
stresses  well  below  the  upper  yield  point.16'17  Such  phenomena 
are  well  known  for  metals  where  microyield  also  occurs  well  below 
the  macroscopic  yield  point18  with  sufficient  mobility  and 
multiplication  rates  for  dislocations  as  can  occur  at  high 
temperatures.  The  microyield  regime  A19  merges  smoothly  into 
stage  II  deformation  behavior. 

The  stress-strain  curves,  with  the  exception  of  the  In-doped 
7 00°C  case,  are  characterized  by  a  stage  II  linear  work 
hardening  regime  of  flow,  followed  by  stage  III  characterized  by  a 
decrease  in  work  hardening  rate,  and  the  beginning  of  the  stage 
IV  linear  hardening  regime  characteristic  of  diamond  cubic 
materials. 14  » 15 , 2  0 

The  results  indicate  hardening  effects  by  indium  in  all 
stages  of  deformation.  The  stage  A  slopes  increase  with 
decreasing  temperature  and  are  larger  at  all  temperatures  for  the 
In-doped  crystals.  The  onset  of  stage  III  occurs  at  larger 
stresses  in  the  In-doped  case,  and  the  slope  in  stage  IV  is 


greater  for  the  In-doped  case.  A  yield  point,  indicative  of 
insufficient  dislocation  mobility/multiplication,  occurs  for  the 
In-doped  crystals  at  700°C  but  not  for  the  undoped  crystals. 

Also,  for  the  700 °C  case,  the  microyield  flow  region  A,  below  the 
yield  point,  has  a  larger  slope  (larger  work  hardening)  than  the 
post-yield,  stage  II,  linear  work  hardening  regime. 

Values  for  the  0.2%  offset  yield  strength  from  the  present 
work  and  that  of  others  are  presented  in  Table  I  and  plotted  in 
Figure  3.  The  data  are  all  converted  to  critical  resolved  shear 
stresses  for  later  comparison  with  single  slip  crystals.  The 
agreement  shown  in  Figure  3  is  seen  to  be  quite  good. 

(B)  [123]  Orientation 

The  specimens  tested  in  [123]  orientation  showed  sharp 
yield  points.  Flow  stress-strain  curves  for  undoped  and  In-doped 
specimens  at  700°C  are  shown  in  Figure  4.  The  single  slip 
orientation  cyrstals  showed  extensive  stage  I  deformation  with  a 
small  yield  drop  at  700°C  (Figure  4),  but  not  at  900  or  1100°C. 

The  pre-yield  slopes  in  the  microyield  region  A  varied  just  as  in 
the  [001]  case,  with  smaller  slopes  for  the  undoped  case  at  each 
temperature  and  with  a  decrease  in  slope  with  increasing 
temperature. 

In  comparison,  the  suppression  of  yield  point  phenomena  at 
temperatures  above  700°C  for  [001]  crystals  (present  work,  ref. 

11)  indicates  a  strong  impediment  to  dislocation  intersection  in 
these  crystals.  This  impediment,  favored  by  the  relatively  low 
stacking  fault  energy  of  -45  mJ/m2  in  these  materials,21 
suppresses  easy  glide  and  contributes  to  the  large  stage  II  slope 


for  the  [001]  crystal.  Neverthless,  the  dislocation  mobility  in 
stage  II  is  sufficient  to  prevent  macroscopic  yield  point 
behavior. 

The  values  for  the  upper  yield  point  as  a  function  of 
temperature,  together  with  results  from  other  work,  are  listed  in 
Table  I  and  plotted  in  Figure  5.  Lower  temperature  data  for 
[100]  specimens  from  the  work  of  Hobgood  et  al.6  and 
Swaminathan  and  Copley9  are  included  because  their  specimens 
also  exhibited  sharp  yield  points  with  little  yield  drop.  This 
[001]  yield  point  behavior  below  700°C  is  associated  with  the 
reduction  in  dislocation  multiplication/mobility  by  the  presence 
of  the  peierls  barrier  and  the  attendant  mechanism  of  double-kink 
dislocation  motion. 

The  present  results  align  well  with  the  other  work  on  the 
same  material.6,9'10  All  results  are  in  agreement  with  regard 
to  the  general  shape  of  the  curves  of  flow  stress  versus 
temperature.  However  the  present  results  are  consistently  higher 
in  stress,  by  about  a  factor  of  two,  than  the  other  extensive 
work  on  single  slip  crystals.10  A  possible  cause  for  this 
difference,  as  discussed  previously®  is  the  presence  of 
-5  x  1017cm-3  of  boron  in  the  present  samples,  a  solute  that  would 
also  be  expected  to  produce  strong  solution  hardening. 

DISCUSSION 

(A)  Compression  Tests 

The  overall  form  of  the  stress-strain  curves  with  stages  II 
to  IV  for  [001]  crystals,  and  stage  I,  as  well,  for  [123]  crystals, 
is  consistent  with  tne  behavior  at  elevated  temperatures  of  the 


much  more  extensively  studied  Si  and  Ge  crystals.14,15,20  The 
analogy  suggests  that  glide  processes  are  important  in 
determining  the  flow  stress  in  stages  I  and  II  while  recovery 
processes  are  important,  in  addition,  in  stages  III  and  IV.  The 
relative  temperature  independence  of  the  critical  resolved  shear 
stress  above  700°C  (Figures  3  and  5),  implies  plateau-type 
behavior**  characteristic  of  athermal  solid  solution  hardening. 
Hence,  the  hardening  indicated  at  yield  and  in  stage  I  is 
consistent  with  the  suggestion7  that  In  should  provide  strong 
solid  solution  hardening  in  Ga1_xlnxAs. 

Below  ~700°C,  the  sharp  rise  in  flow  stress  with  decreasing 
temperature  is  related  to  a  change  in  the  mechanism  of  glide  to 
one  controlled  by  double-kink  nucleation  and  propagation,  again 
analogous  to  Si  and  Ge.  The  hardening  effect  in  this  regime 
could  also  be  related  to  the  elastic  field  of  an  InAs4  unit,7 
in  this  case  effective  through  an  impediment  to  the  propagation 
of  kinks  along  a  dislocation  line.  Because  of  the  apparent 
change  of  mechanism  below  about  700°C,  results  obtained  below 
this  temperature,  e.g.  on  dislocation  velocity,23,24  may  not  be 
applicable  to  behavior  above  700°C. 

The  larger  stress  for  the  onset  of  stage  III  and  the  larger 
slope  in  stage  IV  for  In-doped  crystals  imply  a  decrease  in 
recovery  rate  for  the  In-doped  crystals.  This  could  be  caused, 
in  principle,  by  suppression  of  cross-slip  or  by  retardation  of 
climb.  Several  factors  indicate  that  the  latter  climb  effect  is 
predominant.  First,  transmission  electron  microscopy  studies21 
indicate  that  In-doping,  at  the  level  studied  here,  does  not 


affect  the  stacking  fault  energy,  implying  that  there  should  be 
little  influence  on  the  cross-slip  probability.  Second,  our 
transmission  electron  microscopy  studies25  of  [001]  specimens 
strained  into  stage  IV  show  well  developed  dislocation  networks 
similar  to  structures  observed  for  stage  IV  deformation  of 
Si,  a  case  where  climb  is  known  to  be  the  dominant  recovery 
mechanism.  Third,  stress  relaxation  measurements  in  the  present 
work26  show  consistently  lower  relaxation  rates  for  the  In-doped 
case  at  a  given  stress  over  the  entire  range  of  relaxation, 
implying  a  lower  climb  rate. 

(B)  Comparison  with  Hardness  Results 

Resolved  shear  stress  values  calculated  earlier  by  us  from 
values  of  the  hardness  H  8  are  much  higher  than  the  measured 
flow  stress  values.  Comparative  data  using  the  earlier  estimate 
of  H/6  together  with  4%  offset  flow  stresses  are  listed  in  Table 
I.  The  reason  for  the  discrepancy  is  our  earlier  use  of  an 
empirical  correlation  of  hardness  and  yield  strength  for 
engineering  materials  with  characteristically  low  work  hardening 
rates.  If,  instead,  the  strain  hardening  rates  are  taken  into 
account  using  the  empirical  relationship  developed  by  Cohoo.n  et 
al.,  (H/  ^7,4)  =  3  x  10n  where  n  is  the  strain  hardening 
exponent,  and  °~i is  the  flow  stress  at  7.4%  strain,  the  values 
of  flow  stress  calculated  from  hardness  values  approach  the  flow 
stress  at  4%  strain,  but  still  appreciably  exceed  the  0.2%  offset 
yield  stress.  The  trends  in  hardness  with  temperature  and  the 
hardening  increment  provided  by  indium  as  revealed  in  the  earlier 
work,  however,  do  give  an  indication  of  the  strengthening  effect 
of  In  consistent  with  the  present  work. 
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(C)  Crystal  Growth 

Recent  calculations  by  Jordan  et  al.28  indicate  that 
dislocations  can  be  completely  removed  in  75  mm  diameter  GaAs 
crystals  grown  under  low  thermal  gradient  conditions  if  the  CRSS 
for  plastic  flow  initiation  is  a  factor  of  four  larger  than  a 
GaAs  base  level  value  of  0.6  MPa,  determined  by  extrapolation  to 
the  melting  point.  The  basis  for  the  extrapolation  is  the  data 
for  the  yield  point  of  undoped  GaAs  (ref.  9,  Figure  3).  The 
present  results,  Figure  3  and  5,  suggest  a  value  of  3.3  MPa  for 
the  yield  point  extrapolated  to  the  melting  point  for  the  present 
case  of  In  at  a  level  of  1-2  x  1020cm“2.  Hence,  such  a  doping 
level  of  indium  should  prevent  the  profuse  dislocation  generation 
characteristic  of  stage  I  deformation  according  to  the  theory  of 
Jordon  et  al.28 

However,  the  present  work  and  that  of  others6'10,11  show 
that  microyield  occurs  well  below  the  yield  point  in  compression 
tests  above  700°C  and  particularly  at  1100°C.  These  results, 
together  with  the  in  situ,  X-ray  topographic  studies  in 
tension16'17  at  temperatures  less  than  700  °C  ,  indicate  that 
dislocations,  once  formed,  are  highly  mobile  at  high  temperatures. 
This  suggests  that,  if  GaAs  crystals  are  grown  in  a  dislocation- 
free  condition,  a  major  impediment  to  the  presence  of  dislocations 
must  be  the  nucleation  process.  In  all  of  the  compression  tests 
[present  work,  refs.  6,  9,  10,  11]  the  crystal  end  in  contact 
with  the  loading  platens  is  an  easy  site  for  dislocation 
nucleation  because  of  stress  intensification  at  contact 
asperities  and  because  of  the  weak  singularity  in  the  stress 


field  at  the  edge  of  the  contact  surface  of  the  specimen.  Hence, 
such  data,  while  indicative  of  the  yield  point,  may  not 
correspond  in  microyield  behavior  to  the  crystal  growth  case 
where  nucleation  of  dislocations  would  be  more  difficult.  There 
is  some  indication  of  such  a  possibility  in  the  low  temperature 
tensile  test  observations  of  Tohno  et  al., 16,17  a  case  where  such 
end  effects  are  absent  and  where  dislocation  nucleation  does  not 
occur  at  stresses  of  at  least  twenty  percent  of  the  yield 
stress . 
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SUMMARY 


(A)  Addition  of  In  at  levels  of  1-2  x  1020  atoms/cm3 
results  in  CRSS  values  for  flow  that  are  about  twice  as  large  as 
the  undoped  crystals  throughout  the  stress-strain  curve  and  in 
the  temperature  range  of  700  to  1100°C. 

(B)  The  CRSS  for  a  yield  point  or  for  the  0.2%  offset  yield 
strength  was  weakly  dependent  on  temperature  in  the  temperature 
range  investigated  as  expected  for  a  model  of  athermal  solid 
solution  hardening. 

(C)  The  strength  of  In-doped  crystal  is  sufficient  to 
eliminate  profuse  dislocation  formation  in  a  75  mm  diameter 
crystal  on  the  basis  of  current  theories  for  the  magnitude  of  the 
thermal  stress  during  growth. 

(D)  The  onset  of  stage  III  occurs  at  much  higher  stress 
levels  in  the  In-doped  alloy.  Together  with  observations  of 
dislocation  structures,  this  suggests  that  the  process  of 
dislocation  climb  is  slowed  appreciably  by  In-doping. 
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FIGURE  CAPTIONS 


Figure  1.  Schematic  of  the  compression  testing  jig. 

Figure  2.  Engineering  stress  strain  curves  at  different 
temperatures  for  a)  In-doped  GaAs  and  b)  undoped  GaAs  specimens 
tested  in  the  [001]  orientation. 

Figure  3.  Comparison  of  critical  resolved  shear  stress  values  of 
undoped  and  In-doped  GaAs  from  present  and  other  work  in  the 
[001]  orientation  as  a  function  of  temperature.  Open  squares, 
In-doped  (this  study);  solid  squares,  undoped  (this  study); 
open  triangle,  In-doped  (Ref. 6);  solid  triangle,  undoped  (Ref. 

6);  open  circles,  In-doped  (Ref.  11);  solid  circles,  undoped 
(Ref.  11). 

Figure  4.  Engineering  stress  strain  curves  for  undoped  and  In- 
doped  GaAs  in  the  [123]  orientation. 

Figure  5.  Comparison  of  critical  resolved  shear  stress  values  of 
undoped  and  In-doped  GaAs  from  present  and  other  work  as  a 
function  of  temperature.  Inverted  open  triangle,  In-doped  (this 
study);  inverted  solid  triangle,  undoped  (this  study);  open 
triangle,  In-doped  (Ref.  6);  solid  triangles,  undoped  (Ref.  6)  ; 
open  squares.  In-doped  (Ref.  10);  solid  squares,  undoped  (Ref. 
10);  solid  circles,  undoped  (Ref.  9). 
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Table  I.  Resolved  shear  stress  values  for  ur.dcpec  and  In-d:ped  GaAs. 

Oy  is  upper  yield  point,  oQ  ^  and  oQ  are  offset  yield  strengths 
at  strains  of  0.002  and  0.04,  respectively,  and  H/6  is  an  estimate  fro:  data  for 
hardness  H.  Numbers  in  parentheses  are  references  to  other  work. 
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Figure  1.  Schematic  of  the  compression  testing  jig 


circles,  undoped  (Ref.  11). 
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Figure  4.  Engineering  stress  strain  curves  for  undoped  and  In 
doped  GaAs  in  the  [123]  orientation. 
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TRANSMISSION  ELECTRON  MICROSCOPY  STUDIES  IN  DEFORMED  GaxIni_xAs 
R.  S.  Rai,9  K.  T.  Faber,9  S.  Guruswamy  9  and  J.  P.  Hirth" 

9  Department  of  Ceramic  Engineering,  and  "Department  of  Metallurgical 
Engineering,  The  Ohio  State  University,  Columbus,  Ohio  43210 

The  control  of  dislocation  density  during  the  growth  of  GaAs  and  related 
compounds  is  highly  desirable  for  obtaining  improved  performance  and 
reliability  of  opto-electronic  devices.  Doping  of  single  crystal  GaAs  grown 
by  the  LEC  process  with  Indium  has  been  known  to  reduce  the  dislocation 
density  significantly.  Substitutional  solid-solution  strengthening  of  GaAs 
as  an  InAs  unit  has  been  suggested  to  be  responsible  for  reduction  of 
dislocation  ^density.  3  To  understand  the  mechanism  involved  in  dislocation 
density  reduction,  deformation  tests  have  been  performed  on  [Ogl]  oriented 
Ga  In.  As  single  crystals  in  the  temperature  range  700-1100  C  and  this 
pa$er  reports  some  results  of  the  TIM  characterization  of  dislocations  in 
these  deformed  single  crystals. 

Specimens  of  GaAs  and  Gan  —In.  mAs  for  TEM  observations  were  made  by 
cutting  thin  slices  in  a  disifed  orientation,  followed  by  mechanical 
polishing  and  ion-milling.  These  specimens  were  examined  in  a  JE'U.  JEM  20GCX 
transmission  electron  microscope  at  200KeV  equipped  with  a  double-tilt 
goniometer  stage. 

In  GaAs,  owing  to  its  polarity,  two  types  of  dislocations  with  opposite 
Burgers  vectors ,  termed  0  and  S  ,  can  be  distinguished .  Both  dislocation 
types  were  four.,  in  the  present  study.  Edge  character  dominated  in  the  higher 
temperature  deformed  specimens  in  contrast  with  the  low  temperature  unaoped 
deformed  specimens  where  screw  dislocations  dominate.**  In  the  7G0  deformed 
specimens,  In-doped  and  undoped  both  show  the  presence  of  small  dislocation 
loops , complex  dislocations  networks  and  some  isolated  dislocations.  Only  in 
case  of  the  In-doped  specimens,  straight  dislocations  lying  along  <110> 
directions  are  seen  in  some  regions.  An  example  of  the  dislocation  structure 
in  the  700°C  undoped  specimen  is  shown  in  Fig.1.  In  a  few  places, 
dislocations  are  dissociated  into  two  Shockley  partials.  The  dislocations 
have  been  characterized  using  the  standard  g.b  criterion.  For  partial 
dislocations, the  analysis  was  carried  out  with  the  help  of  weak  beam  pictures 
taken  in  <220>  type  reflections. 

Specimens  deformed  at  1000  and  1100°C  show  quite  relaxed  dislocation 
structures  compareu  to  the  700  C  case.  Hexagonal  networks  of  dislocations, 
parallel  sets  of  dislocations  intersected  by  dislocations  of  different  slip 
systems,  as  well  as  isolated  dislocations  were  found.  Fig. 2  shows  hexagonal 
networks  of  dislocations  in  the  1000  C  deformed  In-doped  specimen. 
Dislocations  of  one  set  are  clearly  dissociated  into  two  Shockley  partials. 
One  interesting  examole  is  shown  in  Fig. 3  where  dislocation  lines  lie  nearly 
along  <11 0>  and  <11 2>  directions.  Most  of  the  dislocations  in  these 
conf igurations  were  found  to  be__  dissociated.  The  dissociation  scheme 

a/2 [ 1 01 i  _ >  a/6  [211]  +  a/6  [112]for  dislocations  A  (marked  in  Fig. 3)  is 

the  expected  one  that  forms  an  intrinsic  stacking  fault. 

From  this  study,  dislocation  structures  in  these  Ga  In^  xAs  alloys  consist  of 
dislocation  loops,  networks  of  dislocations  and  isolated  dislocation  dipoles. 
A  detailed  TEM  study  of  deformed  specimens  is  still  in  progress  and  will  be 
reported  elsewhere.' 
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FIG.  1. — Bright  field  image  ,  g[25o],  foil  orientation  [001]. 

FIG.  2. — Bright  field  image  showing  hexagonal  network  of  dislocations, 
g[022],  foil  orientation  [Oil].  _ 

FIG.  3. — Bright  field  image,  g [ 202 ]  and  foil  orientation  [111]. 
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Appendix  C 


1  Damage  of  coherent  multilayer  structures  by  injection  of  dislocations 

or  cracks 

•  J.  P.  Hirth 

1  Metallurgical  Engineering  Department.  Ohio  State  University.  Columbus.  Ohio  43210 

A.  G.  Evans 

|  Materials  Program.  College  of  Engineering,  University  of  California.  Santa  Barbara.  California  93106 

k  (Received  24  February  1986;  accepted  for  publication  13  June  1986) 

k  Experimental  evidence  indicates  that  coherent  multilayer  structures  (strained  superlattices) 

k  can  be  grown  in  a  damage-free  state  when  the  thickness  of  the  layer  is  less  than  a  critical  value. 

The  resistance  of  such  structures  to  subsequent  damage  by  dislocation  or  crack  injection  is 
examined  in  the  present  study.  The  structures  are  determined  to  exhibit  damage  resistance  that 
I  decreases  as  either  the  layer  thickness  or  the  coherency  strains  increase. 
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I.  INTRODUCTION 

There  has  been  burgeoning  recent  interest1’  in  struc¬ 
tures  composed  of  alternating  layers  of  different  crystals.  Of 
particular  import  are  coherent,  dislocation-free,  multilayer 
structures  (also  called  “strained  superlattices’'),  which  ex¬ 
hibit  interesting  electronic  properties  because  of  the  reduced 
symmetry  associated  with  the  large  coherency  strains.  How¬ 
ever,  multilayer  structures  can  only  be  grown  in  a  disloca¬ 
tion-free  condition  when  the  layers  have  thicknesses  less 
than  acritical  value.5'5  Furthermore,  such  structures  tend  to 
be  unstable  in  the  presence  of  mechanical  and  thermally  in¬ 
duced  stresses,  particularly  stresses  resulting  from  laser  irra¬ 
diation. '’•7  In  the  present  article  the  stability  of  coherent  mul¬ 
tilayer  structures  is  analyzed  with  respect  to  injection  of 
either  dislocations  or  cracks.  For  this  purpose,  the  residual 
coherency-induced  stresses  that  exist  within  a  multilayer 
system  are  first  described.  Aspects  of  dislocation  injection  in 
the  Ga-As-P  system  are  then  evaluated.  Finally,  conditions 
for  the  growth  of  interfacial  edge  cracks  and  the  decohesion 
of  the  extreme  layers  are  estimated. 


with/i  the  shear  modulus  and  v  Poisson's  ratio;  while  for  the 
anisotropic  case 

c  =  cM  +  c,2  —  2(c{2/cn), 

where  c,;  are  the  elastic  constants  referred  to  the  cube  sym 
metry  axes. 

With  reference  to  layers  A  and  B  ( Fig.  1 ),  mechanical 
equilibrium  dictates  that  the  net  force  on  the  system  must  be 
zero,  whereupon, 

aAhA  +  <JBhB  =  0,  .  (2) 

where  h  is  the  layer  thickness.  The  misfit  strains  must  also 
partition  such  that 

e0  =  £*  +  G,=iW<a>,  (3) 

where  a  is  the  lattice  parameter,  Aa  =  aB  —  aA .  and  (a)  su 
(since,  in  all  cases  of  interest.  e„<0.01).  Combining  Eqs. 
( 1 ),  (2),  and  (3)  the  strain  at  equilibrium  is  determined  to 
partition  according  to 

eA  =e0 hBcB/(hAcA  +h„cB ), 

tB  =eJ>.*cA/(hAcA  +hBcB),  (4) 


II.  RESIDUAL  STRESSES 

Coherency  between  layers  in  a  multilayer  system  results 
in  residual  stresses.  The  stresses  are  spatially  uniform  within 
the  plane  of  the  layer,  except  in  the  vicinity  of  edges.  The 
magnitudes  of  the  uniform  stresses  and  the  effects  of  edges 
are  derived  as  a  basis  for  predicting  the  dislocation  and  crack 
injection  conditions. 


with  the  signs  of  eA  and  eB  depending  on  the  magnitudes  of 
a  A  and  a  B . 

B.  Edge  stresses 

Calculation  of  the  edge  stresses  induced  by  the  coheren¬ 
cy  strain  is  achieved  (for  elastically  isotropic  systems)  by 


A.  Uniform  stresses 

Lattice  mismatch  for  coherent  {100}  interfaces  gives 
rise  to  equal  biaxial  stresses  (cr, ,  =  a22=a)  and  strains 
(e,  i  =  f22  =  e„)  in  the  two  layers.  For  a  free  pair  of  layers, 
bending  stresses  would  also  be  present.  However,  these  are 
absent  in  the  complete  multilayer  structures  considered 
here.  The  boundary  condition  an  =  0  dictates  that 

a  =  cf„,  ( 1 ) 

where  for  isotropic  materials' 


c  =  2/i(l  +  v)/(  1  —  »’), 


FIG  I .  A  pair  of  layers  from  a  multiple  layer  structure. 
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imposing  alternating  edge  tractions  that  exactly  cancel  the 
in-plane  stresses  ( Fig.  2 ) .  Such  tractions  create  nonuniform 
stresses  within  the  elastic  half  space."  Furthermore,  the  sym¬ 
metry  of  the  problem  dictates  that  the  largest  localized  edge 
stresses  are  the  interfacial  shear  stresses,  r1"1 .  ( Note  that  no 
normal  stresses  exist  at  the  central  interface,  at  least  when 
the  elastic  constants  are  homogeneous. )  For  one  pair  of  lay¬ 
ers.  within  a  homogeneous  semi-infinite  medium,  integra¬ 
tion  of  the  solution  for  a  line  force  acting  on  an  elastic  half 
space"  gives 

r>"' (2)  =  lce„h1/n{z1  +  hl).  (5) 

The  analogous  result  for  the  stress  at  the  central  interface 
among  .V  pairs  of  such  planes,  is 

(l  +4^(/,2+z-") 

n(r  +  h~)  \ 


x  V  (  ~  i  )"(2h  +  1 )  \ 

[r  +  (n  +  |)-A  ) 


(6) 


Note  that  the  shear  stress  has  a  maximum  at  the  surface 
given  by 


rm"  =  2ceJ  jt  (7) 

and  becomes  negligibly  small  at  z/h  i  4.  The  peak  edge  shear 
stress  is  thus  similar  in  magnitude  to  the  uniform  normal 
stress  (a  =  ce„). 

The  edge  stresses  at  the  interface  between  the  top  layer 
and  the  substrate  are  more  difficult  to  compute,  because  of 
the  presence  of  the  second  free  surface.  Finite  element  re¬ 
sults'  indicate  that  both  normal  and  shear  stresses  exist  at 
this  interface. 


111.  DISLOCATION  INJECTION 

Dislocation  formation  can  be  motivated  by  both  the  co¬ 
herency  and  applied  stresses.  A  lower  bound  estimate  of  the 
layer  thickness  needed  to  form  dislocations,  in  the  presence 
of  these  stresses,  can  be  obtained  by  equating  the  total  force 


no.  2.  Traction  applied  lo  the  surface  in  ordci  to  evaluate  the  shear  stress 
concentration  axz. 


FIG.  3.  View  of  dislocation  configuration,  (a)  parallel  to  the  interfaces  of  an 
A  layer  and  perpendicular  to  the  arms  E  of  the  dipole,  (b)  parallel  to  the 
interfaces  of  the  A  layer  and  parallel  to  the  arms  E  indicating  the  slanted 
configuration  of  the  dipole. 


on  the  dislocation  (the  Peach-Koehler  force)  to  the  line 
tension.  This  simple  calculation  neglects  the  details  of  the 
dislocation  nucleation  process.  A  more  rigorous  approach 
considers  specific  dislocation  loop  configurations  at  nuclea¬ 
tion.  and  determines  the  critical  nucleation  condition.  Both 
approaches  are  afforded  brief  consideration. 

Dislocation  formation  in  Ga-As-P  multilayer  systems 
has  previously  been  considered  '  to  involve  the  creation  of  a 
dipole  on  the  glide  plane  of  GaAs  ( Fig.  3 ) .  Alternatively,  an 
interface  dislocation  may  be  created  at  the  lateral  surface 
(Fig.  4)  Both  configurations  are  examined  in  this  article. 
The  dipole  configuration  is  used  to  illustrate  the  simple  force 
balance  calculation  leading  to  the  prediction  of  a  critical  lay¬ 
er  thickness,  while  the  interface  dislocation  is  used  to  ad¬ 
dress  the  pioblem  of  dislocation  nucleation  at  the  lateral  sur¬ 
face. 


FIG.  4.  Nucleation  of  a  loop  on  the  (001 )  interface. 
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A.  The  Peacn-Koehler  force 

The  force  on  a  dislocation  arising  from  the  coherency 
stresses  in  the  GaAs-GaAs,,  <  P„  5  system  is  conveniently  de¬ 
termined  in  coordinates  fixed  on  the  [001  ]  ( 1  ll)  glide  sys¬ 
tem  in  GaAs  (i'|| [Oil],  j'|| [111],  k'|| [ 2 1 1 ] ).  In  these  co¬ 
ordinates  (Fig.  3),  the  Burgers  vector  is  b  =  b. 0,0,  the  sense 
vector  |  is  parallel  to  k'  and  the  resolved  shear  stress 
is(v6/6)cr.  The  force  on  the  dislocation  is  then  the  product 
of  the  force  pier  unit  length  and  the  segment  length  (>/6/2)/r, 
giving 

Fpk  =ceJ>h/2  (8) 

in  agreement  with  earlier  work.'  When  a  stress  having  a 
shear  component  r  along  with  the  slip  plane  is  subsequently 
applied,  the  force  increases,  such  that 

F=Fpk  +Fa 

=  (bh)(ce^2  +  v/6r/2).  (9) 
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B.  Line  tension  force 

For  the  configuration  of  Fig.  3  the  line  tension  force  on 
the  movable  segment  D  is 


Fs  =  2(/r/47r)ln(v(>Aa/2£),  ( 10) 

where  k  is  the  energy  factor. 1,1  a  is  a  core  cutoff  parameter." 
and  the  factor  2  appears  because  two  fixed  segments  E  act  on 
the  moving  segment  D.  In  the  isotropic  case 

k=  [nb2/(\  —  v)](l  —  vcos2/?),  (11) 


where/?  is  the  angle  between  b  and  |.  In  the  anisotropic  case, 
because  the  Burgers  vector  has  a  screw  component 
by  =  j  [  1 10]  and  two  edge  components  b ,  =([001]  and 
f»,  =  ][  1 10],  can  be  determined  as 

k  =  (K,b]  +  K:b\  +  Kyb\),  (12) 

where 


«,  =  (?[,  +c 


»)(  (?il  ~C|*)C“  ), 
\c;2  (c;,  +c,2  +  2c 

44)  / 


Ki  —  (^22/^11  )^*|» 

K,=  (C^C^)1'2, 

with  c'22  =cM  +H/2,  c[,  =  (c,,cj,  )l/2,  c^=c44 

—  H  / 2,  and  H  =  2c44  +  c,;  —  cn. 


C.  The  critical  thickness 

The  critical  layer  thickness  hc,  above  which  dislocation 
dipole  generation  becomes  feasible  in  the  absence  of  an  ap¬ 
plied  stress,  is  determined  by  the  condition  Fpt  =  F, .  For 
Ga-As-P,  values  of  this  thickness  are  obtained  using  the 
following  quantities:  a  =  3,  typical  of  compound  semicon¬ 
ductors,"  Voigt  average  values  for  /z  =  48  GPa  and 
v  =  0.23,"  A  a  =  0.020  nm  and  a  —  0.565  nm  for  GaAs. 14 
The  results  (Table  I)  are  in  reasonable  agreement  with  ex¬ 
periment,  indicating  that  dipole  formation  in  GaAs  does  not 
appear  to  be  severely  nucleation  limited.  Moreover,  the  iso¬ 
tropic  approximation  gives  results  that  do  not  differ  drasti¬ 
cally  from  those  for  the  anisotropic  approximation.  It  is  also 


TABLE  I.  Calculated  values  of c,  k.  and  h , 


Case 

c  (GPa) 

K 

h,  ( nm ) 

isotropic 

153 

58.5 

30 

anisotropic 

124 

57.7 

39 

expenment 

*** 

35 

evident  from  Eq.  (9)  that,  in  the  presence  of  an  applied 
stress,  the  critical  thickness  is  reduced  to 

h*  =  hc  [  1  —  x/6(r/ct0)  J.  (13) 


D.  Dislocation  nucleation 

Dislocation  nucleation  occurs  most  readily  at  the  free 
surfaces,  where  image  effects  reduce  the  free  energy  of  for¬ 
mation  of  a  critical  nucleus.15  Moreover,  only  a  small  error 
( less  than  5% )  is  introduced  by  approximating  the  loop  as  a 
semicircle  of  radius  r.  For  an  interface  dislocation 
=  i [  1  To ] .  forming  on  the  (001)  interface,  the  resolved 
shear  stress  arising  from  applied  forces  performs  work , 

W‘=-wrb,T/  2.  (14) 

The  corresponding  elastic  energy  increase  in  the  presence  ot 
a  semicircular  dislocation  is15 

U  =  [  (2  —  v)^b  )r/ 4(  1  —  v)  ]  [Inf  8 ra/b, )  —  2]  =Sr, ( 15 

where  S  is  an  effective  line  tension.  However,  in  the  presence 
of  coherency  strains,  there  is  also  an  elastic  energy  release 
The  strain  energy  removed  per  dislocation  is  given  by  the 
product  of  the  strain  energy  density  and  the  stress  relief  vol¬ 
ume.  For  a  semicircular  loop  forming  at  the  surface,  the 
coherency  contribution  to  the  nucleus  formation  energy  W 
can  be  evaluated  from  the  edge  stress  [Eq  (7)  ]  and  the  loop 
area  trr/ 2,  by  recognizing  that  the  zone  of  alternating  strain 
has  a  characteristic  dimension  h.  Then. 

Wc  ~  (  —  cjr5f£/ir:/64)=  —  \ftrrh /2  (16) 

where  M  =  trcei/ 32.  The  total  free  energy  accompanying 
loop  formation  is 

A  G=W°  +  W'  +U.  (17) 

Hence,  by  setting  dAG  /dr  —  0  and  applying  standard  nu¬ 
cleation  theory,15  the  critical  stress  for  dislocation  nuclea¬ 
tion  r*  becomes 

r*=(SV406,A-7')  -  (Mh/b,) 

=(S2/Wb,kT)[  \  —  h  /h*],  (18) 

where  h  *=S2/ 40  MkT  is  the  film  thickness  at  which  the 
barrier  to  nucleation  vanishes.  Typical  values  of  r*  for  un¬ 
strained  crystals15  (0.01-0.05  //m),  suggest  that  h  *  ranges 
between  7  and  40 //m.  Thus,  h  *  appreciably  exceeds  hc.  indi¬ 
cating  that  nucleation  can  represent  a  substantial  restriction 
on  dislocation  injection.  Dislocation  formation  at  h  thus 
relies  upon  the  pre-existence  of  either  threading  dislocations 
or  of  edge  loops.  Furthermore,  it  is  evident  that  the  vulnera¬ 
bility  to  dislocation  injection  increases  as  either  the  misfit 
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FIG  6.  Trends  in  siress  intensity  factor  with  crack  length  lor  subsurface 
edge  cracks. 


tl 


strain,  the  layer  thickness,  or  the  imposed  stress  increases. 

IV.  CRACK  FORMATION 

The  crack  injection  problem  to  be  considered  involves 
the  formation  of  edge  cracks  between  layers,  caused  by  the 
concentrated  coherency  strains  that  develop  in  the  presence 
of  edges.  The  problem  has  some  analogy  with  delamination 
cracks  that  initiate  from  edges  and  holes  in  laminated  fiber 
composites. 16  The  nature  of  the  problem  is  illustrated  by 
considering  a  system  in  which  the  layers  have  equal  thick¬ 
ness  and  have  the  same,  isotropic  elastic  constants. 

A.  The  crack  driving  force 

For  elastically  isotrupic  systems,  the  residual  ei'ze 
stresses  can  be  used  to  compute  stress  intensity  factors  K  tor 
various  crack  configurations.  Two  bounding  problems  are 
presented  to  illustrate  potential  modes  of  cracking:  interfa- 
ciai  edge  cracks  near  the  center  of  an  infinite  set  of  layers 
(Fig.  5)  and  a  subsurface  edge  crack  at  the  first  interface 
l  Fig.  6).  For  the  centra/  edge  crack,  a  mode  II  stress  intensi¬ 
ty.  K.  exists,  which  can  be  determined  from  the  stress  field 
by  use  of  a  superposition  procedure.’’ 

...  ,i9) 

\a  J"  (z~  +  h')\  1  —  (z/a)- 

wliere  F(z'a)  =  [1  -  (z/a)']  [0.29  -  0.39(z/a):  +  0.77 
tz/a)*  —  0.99(*/a)6  +  0.51  (i/a)"]  and  a  is  the  crack 
length.  Integration  ofEq.  (19)  gives 

A',,  ~  (2.4ce„vA  /v  ir){(a/h)/[  1  +  (a/h)z]}'n.  (20) 

Trends  in  K  „  /ce,t\h  with  normalized  crack  length  a/h  are 
plotted  in  Fig.  5.  A  maximum  exists  at  a/h  =  1  given  by 

kuzz(2A/\2ir)ce„\7i.  (21) 

Note  that  K ,,  — 0  as  a/h  becomes  large  (  —  -r.  ). 

For  the  subsurface  edge  crack,  numerical  procedures  are 
needed  to  determine  the  driving  force  for  cracks  having  im¬ 
mediate  size.  0.2Sa//tS5.  However,  for  large  a/h.  the 
crack  driving  force  may  be  estimated  from  the  strain  energy 


FIG.  5  Trends  in  normalized  stress  intensity  factors  with  crack  length  for 
central  edge  cracks. 
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release  in  the  upper  layer.  Specifically,  the  strain  energy  in 
the  layer  above  the  crack  is  eliminated  as  the  crack  pro¬ 
gresses  and  thus  constitutes  a  strain  energy  release  rate  Fur¬ 
thermore,  the  strain  energy  concentrated  around  the  crack 
tip  simply  advances  with  the  crack  (i.e..  remains  constant 
during  crack  propagation)  and  thus,  does  not  contribute  to 
the  energy  release  rate  G.  Consequently.  G  can  be  deduced 
from  the  strain  energy  U  in  the  layer  as 

G  =  2  —  -■  -  ( 1  —  v)ha  =  (ce„)zh(  1  —  v)/E. 

aa  da  2  E 

(22) 

In  view  of  the  equivalence.  G  =  K  :( 1  —  tf  )/E.  the  effective 
stress  intensity  factor  becomes 

K/Cf0ih  =  1/v  1  4-  v.  (23) 

This  analysis  does  not  reveal  whether  K  is  single  or  mixed 
mode.  Hoyvever,  recent  calculations'"  have  revealed  that 
both  mode  I  and  mode  II  are  involved,  such  that 

Ktl/K,  =  1.4.  (24) 

Evidently,  Eq.  (23)  reveals  that  K  is  independent  of 
crack  length  for  large  cracks.  Furthermore,  numerical  anal¬ 
ysis'*  reveals  that  K  asymptotically  approaches  the  long 
crack  limit,  as  depicted  in  Fig.  6. 


B.  Crack  growth 

Assessment  of  the  incidence  of  edge  cracking  requires 
that  the  crack  driving  force  be  compared  with  the  crack 
growth  resistance  of  the  interface.  K (typically  in  the 
range  of  0.5-1  MPaytn  for  semiconductors).1'' However,  as 
yet,  a  generalized  criterion  for  mixed  mode  interfacial  frac¬ 
ture  has  not  been  determined.  The  simplest,  reasonable  pos¬ 
tulate  is  that  interfacial  fracture  occurs  when  K  =  A'"".  By- 
invoking  Eqs.  (21 )  and  (23),  one  then  sees  that  fracture  may 
be  completely  averted  (i.e.,  K  never  exceeds  A';nl )  when. 

ceJhSKT-  (25) 

Absolute  assurance  that  the  coherency  strain  does  not  cause 
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FIG.  7.  The  effect  of  applied  loads  r„  on  the  stress  intensity  factor  for  inter- 
facial  edge  cracks. 


cracking  thus  requires  that  the  layer  thickness  should  not 
exceed  a  maximum  value 

h*  =(K™/ce„)'-.  1 26) 

For  the  properties  of  Ga-As-P  specified  above  and  using 
A'"":;  0.5  MPa/m,  h  *  is  predicted  as  40  nm.  similar  in  mag¬ 
nitude  to  h,  (Table  I).  However,  the  prediction  is  sensitive 
to  variations  in  A' ml.  In  order  to  determine  the  specific  value 
i  ."layer  thickness  ( above  h  * )  at  which  cracks  actually  form, 
independent  knowledge  of  the  magnitude  of  interfacial  edge 
defects  would  be  required.  Furthermore,  whereas  central 
edge  cracks  would  propagate  inward  by  at  most  a  few  times 
the  layer  thickness.  ( because  of  the  diminishing  AT ,,  ala  >  It, 
Fig  6),  subsurface  edge  cracks  could  cause  complete  separa¬ 
tion  of  the  top  layer  (Fig.  6). 

When  an  additional  stress  is  imposed  on  the  structure, 
the  cracking  tendency  would  generally  be  enhanced.  The 
details  would  depend  on  the  stress  distribution  in  the  vicinity 
of  the  edges.  For  example,  for  the  simple  case  of  a  uniform 
edge  shear  stress  rR  the  stress  intensity  at  central  edge 
cracks  would  become 


_  /  aV”  (  2.4/N  r  +  vfjj 

ce„sh  '  h  J  \  [  I  +  (a/h)2 ] 1  :  ce„ 


(27) 


The  normalised  stress  intensity  is  thus  modified  by  the 
applied  stress,  as  depicted  in  Fig.  7.  Of  particular  interest  is 
the  minimum  value  of  A',,,  which  essentially  repiesents  a 
barrier  to  catastrophic  fracture.  Specifically,  when  this  mini¬ 
mum  value  of  A' ,,  exceeds  A' complete  fracture  of  the  sys¬ 
tem  becomes  possible. 


V.  CONCLUDING  REMARKS 

The  results  presented  in  this  article  regai  ding  disloca¬ 
tion  injection  and  edge  cracking  in  coherent  layered  struc¬ 
tures.  provide  preliminary  numerical  estimates  pertinent  to 
the  avoidance  of  deleterious  damage  in  such  structures.  The 
estimates  are  independent  of  the  overall  size  of  the  wafer, 
provided  that  the  wafer  diameter  substantially  exceeds  the 
layer  thicknesses  and  that  many  layers  are  involved.  The 
results  show  that  such  structures  are  more  vulnerable  than 
unstrained  phases,  with  respect  to  both  dislocation  and 
crack  injection.  The  susceptibility  to  damage  increases  with 
increasing  lattice  mismatch  (increasing  coherency  strain) 
and  with  increasing  layer  thickness.  The  present  analysis  of 
dislocation  vulnerability  has  neglected  the  lattice  resistance 
to  defect  motion  ( the  Peierls  stress  or  one  of  several  possible 
pinning  effects ) .  Such  factors  should  be  important  at  room 
temperature  (in  Ga-As-P.  for  example),  resulting  in  en¬ 
hanced  damage  tolerance  at  such  temperatures  However, 
these  contributions  to  dislocation  resistance  decrease  in  ini 
portance  with  increasing  temperature. 

ACKNOWLEDGMENTS 

This  research  was  supported  in  part  ,  J  P.H. )  by  the  Air 
Force  Office  of  Scientific  Research  under  Grant  Nr 
F49620-85-C-0129  and  in  part  (A  G  E.)  by  the  Office  i 
Naval  Research  under  Contract  No.  N00014-85-K-0362 

‘L.  Esaki  and  R.  Tsu.  IBM  J  Res.  Dev  14  686  (  l»70). 

K.  Ploog  and  G.  H.  Dbhler.  Adv  Phys.  32.  285  (  1983). 

*J  W  Matthews  and  A  E.  Blakeslee.  J  Cryst.  Growth  27.  1 18  ( 1974* 

4J  W  Matthews  and  A.  E.  Blakeslee.  J  Crvst  Growth  29.  273  ( 197? ) 

'J  W  Matthews  and  A.  E  Blakeslee.  J  Cryst.  Growth  32.  2o5  :  W7o 

*N.  Holonyak.  Jr..  R  M.  Kolbas.  W  D  Latdig.  M.  Altarelli.  R  D  Dupa: 

and  P  D  Dapkus.  Appl  Phys.  Lett  34.  502  i  1 97*1 . 

'G.  C.  Osbourn.  J.  Appl  Phys.  53.  1586  ( 1982) 

*S  Timoshenko  and  J  N.  Goodier.  Theory'  of  Elasticity  (McGraw-H.! 
New  York.  1950) 

'M.  D.  Drory  and  A.  G.  Evans,  unpublished  results. 

,nJ.  P  Hirth  and  J  Lothe.  Theory  of  Dislocations,  2nd  ed  (Wiley.  New 

York.  1982).  p.  91. 

"J.  P  Hirth  and  J.  Lothe.  Theory  of  Dislocations.  2nd  ed.  'Wiley.  New 
York.  1982).  p.  231. 

,:J.  P  Hirth  and  J  Lothe.  Theory  of  Dislocations.  2nd  ed.  (Wiley.  New 
York.  1982).  p  453. 

MK.  Hellwege  and  A.  M.  Hellwege.  Landoll-Bbrnstein  Tables  (Springer 
Berlin.  1969),  Vol.  II1-2. 

,JR  W  G.  WyckoflT,  Crystal  Structures  (W’iley.  New  York.  1963). 

'“J  P  Hirth.  in  The  Relation  Bern  een  the  Structure  and  Mechanical  Prone • 
ties  of  Metals  (H  M.  Stationary  Office.  London.  1963).  p  2 17. 
lf,R.  B  Pipes  and  N.  J  Pagano.  J  Compos  Mater  4.  204  ( 1970) 

,7G.  C.  Sih.  Handbook  of  Stress  Intensity  Factors  (  Lehigh  University  L« 
high.  PA.  1973). 

'*M.  D.  Thouless,  A  G.  Evans.  M  F  Ashby,  and  J  W.  Hutchinson  un¬ 
published). 

WB.  R.  Lawn  and  T  R.  Wilshaw.  Fracture  of  Brittle  Materials  (Cambridge 
University.  Cambridge  Lngland.  1975) 


s 


s 


% 


■: 

■■ 

n 

C 


ft 

1 

2 


J 

■ 


S 

! 

* 

* 

► 


K 

1 

l"  2376  J.  Appl  Phys.,  Vol.  60.  No  7. 1  October  1 986 


J  P  Hirth  and  A  G.  Evans 


2376 


Miiwsraontx  yaaaeei\x\xiM  /i  vrrt  vi  vwf*iflQnu7\w.x  imaxaaamw  s.s.~  AMI  hflMGflMBHMiMQII MMBG 


